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We calculate the current-voltage characteristic of a homogeneously strained metallic carbon nan- 
otube adsorbed on a substrate. The strain generates a gap in the energy spectrum leading to a 
reduction of the current. In the elastic regime, the current- voltage characteristic shows a large neg- 
ative differential conductance at bias voltages of around > 0.17V. We discuss the implications for 
the current in the superelongated regime. 
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I. INTRODUCTION 

Carbon nanotubes have many electronic and mechan- 
ical applications. Their ability to sustain high currents 
before breaking was observed in a number of experiments 
on metallic single-walled nanotubes (SWNT) adsorbed 
on a substrate [2-3|. Corresponding theoretical work 
was done in Refs. |j, 0. With respect to the mechani- 
cal properties, carbon nanotubes belong to the strongest 
and stiffest materials discovered so far. Under stress, a 
tube will first expand elastically. After reaching the yield 
point, it undergoes plastic deformations, leaving a per- 
manent deformation. This point was measured at strains 
around 5 — 10% for nanotubes with diameter 1 — 3nm 
depending on the chirality [1, 0] ■ 

When the effective temperature of the nanotube caused 
by the intense electron-phonon scattering at a high bias 
voltage (we assume in this paper that the system is not 
externally heated) exceeds the activation barriers of de- 
fects, large spiral-like defect stripes are created along the 
tube Q. These defect stripes correspond to missing lat- 
tice plane pieces in the crystal. The missing atoms of 
these plane pieces are rearranged at the end points of 
the nanotube in such a way that a superelongation of 
SWNTs of up to 270% before cracking is reached. In 
Ref. [t| we have developed a theory for this effect in a 
simplified crystal model. Finally, when the effective tem- 
perature due to the applied bias voltage is lower than 
the activation barrier the system responds to the exter- 
nal stress with large non-homogeneous defect pile ups 
or immediately bond-breaking depending on the physi- 
cal surrounding [l(| El ■ This leads then to necking of 
the SWNT. Here, we will not address the current- voltage 
characteristic in this regime. 

The purpose of this paper is to analyse the electrical 
properties in the stress ranges where nanotubes show a 
homogeneous behavior in space, i.e., in the elastic regime 
of small elongation, and in the superelongation regime. 
We will first discuss the elastic regime, and the supere- 



longation regime will be discussed at the end. 

If subjected to an external homogeneous stress, a nan- 
otube responds with a homogeneous strain. This results 
in modifications in the bond lengths and leads, via the 
electron-phonon interaction, to a gap in the electron spec- 
trum of metallic SWNTs. This strongly modifies the 
current-voltage characteristic, and represents, therefore, 
a promising effect for building sensitive stress sensors. 

On a substrate, the current increases with voltage, 
whereas for suspended SWNTs it decreases, i.e., it shows 
a small negative differential conductance at high voltages 
[13] ■ Similarly, we will find here a negative differential 
conductance which is now very large for a SWNT ad- 
sorbed on a substrate under external stress at bias volt- 
ages of around 0.17V. 

The energy levels of electrons in a nanotube form one- 
dimensional bands in the graphene Brillouin zone around 
the K and K' points. For metallic SWNTs, there are 
two energy bands corresponding to right- and left-moving 
electrons which cross at these points. We assume here 
that the diameter of the nanotube is so small that elec- 
tron excitations to higher bands are negligible. For exam- 
ple this approximation is valid at bias voltages U < 2V 
for nanotubes with diameter 2nm. These diameters are 
typical in present day current-voltage experiments [IHH . 
But we expect the validity of the approximation even be- 
yond this value since electron scattering to higher bands 
is effectively forward due to the band edges where elec- 
trons reverse their direction. Note that we use in our 
numerical calculations below scattering values for nan- 
otubes with a diameter of 2nm |{|. 

In order to calculate the current-voltage characteris- 
tic at high bias voltage we will use the semi-classical 
Boltzmann equation for electrons and phonons. The var- 
ious scattering mechanisms between electrons, phonons 
and impurities will be discussed first in Sect. II. Then 
we will discuss in Sect. Ill the effect of strain on the 
wavefunctions and the energy spectrum of the electrons 
in the nanotube. This leads to modifications in the 
electron-phonon scattering times. Section IV introduces 
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the different Boltzmann equations for the electrons and 
phonons which we use in order to calculate the current- 
voltage characteristic. The results of this calculation will 
be discuss in Sect. V ending up with a conclusion in Sect. 
VI. 



II. SCATTERING PROCESSES 

The method we use here to calculate the current- 
voltage characteristic of metallic nanotubes is based 
on the semi-classical Boltzmann equation. Within this 
method quantum interference corrections to the con- 
ductance are not taken into account It was re- 
cently shown numerically that this correction to the 
conductivity is negligible above room temperature for 
single-walled carbon nanotubes without structural de- 
fects due to phonon scattering decoherence mechanisms 
fl~H . This temperature is immediately reached through 
electron-optical phonon and further optical phonon- 
acoustic phonon scattering in nanotubes under high bias 
voltage. 

At low voltages innervalley elastic scattering, which 
consists of acoustic phonon scattering (in the quasi- 
elastic limit) and impurity scattering, is most relevant. 
We neglect electron-electron collisions. At higher volt- 
ages inelastic scattering between the electrons and optical 
phonons, which then leads to intervalley scattering, be- 
comes relevant. This scattering process leads at high bias 
voltage to the creation of hot optical phonons, which are 
then no longer in thermal equilibrium with the environ- 
ment . This leads us to the conclusion that we have to 
consider besides the electrons also phonons in the semi- 
classical Boltzmann approach. The hopping between an 
electron in the K valley and the other at the K'-valley is 
mediated by zone-boundary optical A[ K-phonons where 
only Kekule type of lattice distortions couple to the elec- 
tronic system |15| . On the other hand jumps of electrons 
within the same valley are mediated by zone-center Ei g 
T-phonons [|[ . The optical phonons decay to lower lying 
secondary acoustic phonons by phonon-phonon scatter- 
ing These then decay to other acoustic phonons or 
phonons in the substrate and to a smaller amount also 
in the leads. We show in Fig. 1 this electron-phonon 
relaxation path. 



III. STRAIN EFFECTS ON THE ELECTRON 
SYSTEM 

By using the electron-phonon interaction Hamiltonian 
for acoustic phonons as in Ref. [IH we have calculated the 
lowest-band cigenfunctions and eigenvalues for a metallic 
nanotube under stress. For electrons in the K valley, 
these are given by (letting the axis of the SWNT point 
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FIG. 1: Main electron-phonon relaxation path of metallic 
SWNTs at high bias voltage where electron-optical phonon 
scattering is dominant. Innervalley electron jumps shown 
as solid doublearrows in the figure are mediated by Ei Y- 
phonons. Intervalley electron jumps mediated by A[ K- 
phonons are shown by dashed doublearrows. 



in the y-direction) 

Kk(y) = e Mk+kyo) (^), (i) 

where s = 1 for the conduction band , s = — 1 for the va- 
lence band and k is the momentum in y-direction. More- 
over, we used the abbreviation k± = k xa ± ik according 
to 

k X(> + ik yo = -^-e l3ri [(u xx -u yy ) + i2u xy ] . (2) 
nvp 

Here vp is the electron velocity, u%j the strain tensor, 
and (72 ~ 2eV [l5| the electron-phonon coupling con- 
stant. The parameter 77 in k XOtVo is the chiral angle of the 
SWNT. For electrons in the K' valley we obtain *f fe (y) 
with the substitution 32 — > —92 and k Xo — ik — > k Xa + ik. 
The energy eigenvalues are 

e(s, k) = shvFyJk% a + k 2 , (3) 

Strain makes k Xa ^ and opens a gap A = 2hvp\k Xo \ in 
the electron spectrum of a formerly metallic SWNT. This 
is consistent with numerical calculations and experiments 
[l7l - l2l| . This spontaneous gap generation is not seen 
by computer calculations in graphene [221] . The energy 
eigenvalues e(s, k) show an important exception: The 
spectrum remains gapless 0, [24| for a metallic SWNT 
with r\ = it/ 6 known as armchair nanotube, where the 
homogeneous axial stress obeys u xy = 0. From this and 
the discussion below it is clear that armchair nanotubes 
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show only a small effect on the current- voltage character- 
istic by an applied external stress. For zigzag SWNTs, 
for which one has rj = 0, we obtain the maximal gap 
value for a given applied axial external stress. 

Let us now discuss the current-voltage characteristic 
due to strain. For a formerly metallic nanotube, we ob- 
tain that the energy band and the electron-phonon scat- 
tering times are modified. Indeed, both depend effec- 
tively on the gap parameter A which is itself a function 
of the strain and the chiral angle of the nanotube. We ne- 
glect here modifications of the eigenfunctions above due 
to k Vo ^ as this leads effectively only to small mo- 
mentum shifts of the zone boundary phoiion velocities. 
The corresponding modifications of the current-voltage 
characteristic are small due to the fact that the effective 
phonon mean free path is much smaller than the SWNT 
length. This has been checked numerically. We now dis- 
cuss the modifications of the effective electron-phonon 
scattering times due to the wavefunction modifications 
under stress. 

We restrict ourselves first to the K valley and zone- 
center phonon scattering. The results for intcrvalley 
scattering are mentioned below Eq. (JUJ). In contrast 
to the current- voltage characteristic of SWNTs without 
band gap [f| we have to take backward scattering (inner 
and intcrband) and interband forward scattering into ac- 
count. Intcrband scattering leads, for example, to the 
result that the conduction band can be filled with elec- 
trons by phonon absorption or depleted by phonon emis- 
sion. Therefore, we should take into account that we have 
from the beginning electrons (holes) injected in the con- 
duction (valence) band in the SWNT by a non-zero lead 
temperature. Without intcrband scattering, only these 
electrons and the corresponding holes contribute to the 
current (see the accompanying discussion to Eq. (fl~3j) ) . 

The electronic zone-center phonon scattering Hamilto- 
nian is given by [25j 
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a \ u y (r) - iu x (r) 



u y {v) + iu x {r) 




(4) 



where a is the nearest-neighbour distance in the lattice 
and g3 was determined in Ref. Q by density functional 
methods for metallic SWNTs. 

Next we calculate the momentum-dependent electron- 
phonon transition matrix elements K^f/J-ffintl^^fc/)) 2 , 
which are proportional to the effective electron-phonon 
scattering times. In the following, we mention for which 
type of phonons, i.e. longitudinal (u y ^ and u x = 0) 
or transversal (u y — and u x ^ 0), the matrix clement 



|(*^|^t|*f 



is non-vanishing: 



l(*f fc >ol^int|^<o)| 2 «|^/|fc + l + fc-/l^ir«4, (5) 

4. (6) 



A" 



^O.fc^Ol^intl^s^O.fc'^O/ 



The matrix element ([5]) comes from interband forward 
scattering, which means s ^ s' , with transverse optical 
phonons {u x ^ 0). For innerband backward scattering, 



i.e. s = s', this non-zero matrix element is mediated 
by longitudinal optical phonons (u y ^ 0). The non-zero 
matrix element ([5]) come from interband backward scat- 
tering with longitudinal optical phonons (u y ^ 0) if the 
gap is small (\k Xa \ <C |fc|, |fc'|), and from transversal opti- 
cal phonons (u x ^ 0) if the gap is large (|/c|, \ k'\ <C \k Xo |). 
Note that ([5]) and ^ become exact for large and small 
gaps and are approximate for intermediate gap values. 
Note also that in the present approximation, (|5|), (J6)), 
the backward scattering matrix elements do not depend 



on the stress, i.e., we can set k x 



0. 



IV. CURRENT- VOLTAGE CHARACTERISTIC 

In section II we introduced a two-valley model for the 
scattering mechanisms of the SWNT under bias voltage. 
In the following, we use the approximation of an effec- 
tive one-valley model coupled to only one sort of optical 
phonons. In fact, this was already used by us in Ref. [l6| 
to calculate analytically the current-voltage characteris- 
tic of unstrained SWNTs with non-reflecting leads. This 
is a simplification of the two- valley Boltzmann approach 
where scattering between the valleys is mediated by 
zone-boundary and zone-center optical phonons [H, Ha | 
which we discussed in Sect. II. The latter model leads 
to current-voltage characteristic of unstrained SWNTs 
0, [IH which are in good accordance with experiments. 
Our one- valley approximation is justified by the fact that 
the optical frequencies of the zone-center phonons and 
zone-boundary phonons, and also the electron-phonon 
coupling in the phononic sectors are rather similar. The 
electron-optical phonon interaction in this model leads 
to an effective inverse scattering time l/r op for the elec- 
trons and 2s p /r ep for the phonons with parameters found 
in Ref. [H 

We denote by /i(fc,cc,<) (fn(k 7 x 7 t)) the left 
(right) moving electron distribution functions, and by 
n : ' s (k,x,t), n bs (k,x,t) the distribution function of the 
phonons mediating forward or backward scattering of 
electrons, respectively. The time evolution of the 
electrons is governed by the semi-classical Boltzmann- 
equation 



eE 

d t T v e (k)d x + —d k 
n 



fh/R= [dtfL/n] 



(7) 



with v e (k) = \de(s,k)/d(hk)\, {d t f L /R\c ~ [dt/z,/ii]e + 
+ [9 t f L / R ]f s + [d t f L / R }bs, and elastic scattering 

[dth\e = (v F /le)[f R {k) - f L (k)] D(e(s, k)) (8) 

consisting of acoustic phonon scattering (in the quasi- 
elastic limit) and impurity scattering where l e is the elas- 
tic scattering mean free path. We denote by v(e(s, k)) 
the density of states of the left or right-moving electrons, 
i.e. v(e(s 1 k)) — l/\de(s, k)/dk\ and v = v(e)KvF is the 
dimcnsionlcss density of states. 
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The time evolution of optical phonons is given by 
[d t + v op (k)d x ] n bs 'f s = \d t n hs 'f s ] + \d t n bs / fs . 

J c L J osc 

_ (9) 

Here v op is the optical phonon velocity [16J. [d t n fcs /^ s ] c 
is a phonon-electron scattering term consisting of the in- 
teraction with backward or forward scattered electrons. 
[d t n bs / fs ] oc is a thermal phonon relaxation term such 



that the coupled electron-phonon system is not heated 
up by applying large voltages on the SWNT. 

Scattering of phonons with electrons leads to the fol- 
lowing scattering contributions (backward scattering and 
forward interband scattering) in the electronic Boltz- 
mann equation ([7} and in the phononic Boltzmann equa- 
tion © (we restrict ourselves to k > phonons in (fTTj) 
for simplicity) 



[dtf L ] bs/fs (k,x) 



(10) 



-L^( c +) r *-//-( c +)|[„W/»( A; + ja :) + l]f R/L {k R/L (e+))[l f L (k L (e))] n bs /f s (k+, x)[l f R/L (k R/L (e+))]f L (k L (e)) 
+ He-)r bs/f %e-)^f s (-k-,x)f R/ ^k R/L (e-))[l-^^ 



d t n bs / fs 



(k,x)=2— 



± 



v{e{kt)) v{e{kl)) 



n bs /f°(k,x) + lU f R (k+)[l f L/R (k- /R )} 



(11) 



+ h/R(-K /R )[l - f R (-k+)} \ - n bs lf s {k,x)\ f L/R {k- L/R )[l - f R (k+)} + f R (-k+)[l - f L/R {-k- L/R )] 

I 



with the abbreviations k^ = fci^e^) — fci(e), 
in pll|) . In (fTT|) we used the abbreviation 



K R/L - ^ 2 



2v F \ 



k 2 



-k 2 



e±hu> 



(12) 



where the right hand side of Eq. (|12| is independent of 
i?, L and fc is the phonon momentum. The numerical con- 
stant s p is given by s p «0.67 The frequency of the 
optical phonon is denoted by w. The factors r'' s /'' ?s (e ± ) 
within our one-valley model are determined by the in- 
fluence of zone-boundary K-phonon scattering on the ef- 
fective electron-phonon scattering time by carrying out a 
similar overlap calculation for these phonons as was done 
for zone-center T phonons below Eq. ([6|). Here we use the 
zone-boundary phonon scattering Hamiltonian in Ref. [26| 
together with similar approximations used below Eq. ([6]). 
This leads to r bs (e ± ) w 0.7(l/i>(e ± ) + l/i>(e)) 2 /4 + 0.3 for 
backward scattering and r^ s w 0.3 for interband forward 
scattering fl6| . 

We restrict ourselves here to short SWNTs with 
lengths between 50nm where phonon scattering becomes 
relevant and 500nm. For the thermal phonon relaxation 
term [dtn bs ^ s ] osc we use an expression which takes into 
account also the second generation acoustic phonons with 
parameters discussed in the case of the metallic SWNTs 
without external stress in Ref. [l6l Since \1/D(e(k^)) — 
l/i>(e(fc^))|~ 1 ^ 1 in (fTTj) we can disregard, in the cal- 
culation of Eq. the term [dtn^ s ] osc , and the term 
v op (k)d x nf s for the forward scattering phonons. To inte- 



grate the Boltzmann equations (|7|)- (fTT]) we use a gener- 
alisation of a discretized time-splitting method described 
in Ref. [ID for the stress-free case. Due to the change in 
the density of states coming from the external stress, the 
momentum and space-grid is no longer equidistant in po- 
sition space but still equidistant in momentum space [27| ■ 
The large prefactor \l/v{e{k+)) - l/£>(e(fc^))|- 1 in CEJJ) 
allows us to determine 9 t n^ s in every iteration process 
by (fTTj) with [dtnl s ~\ c (k , x) = 0. This dynamics leads 
approximately to the same stationary solution of (JTJ) , ((9j) 
as in the case when using the full phonon scattering dy- 
namics (fTTj) . 



V. DISCUSSION 

The upper panel in Fig. 2 with the continuous curves 
shows the current-voltage characteristic calculated with 
the help of ©-(HH) for a SWNT of length L = 300nm 
and different values of the gap energies A. Note that 
the black circles are calculated for gap energies A — > 
just after the gap opening. The recorded data points are 
restricted by the demand that the momentum grid, the 
energy values A and optical phonon frequencies uj are 
commensurate. By comparing the current values / for 
voltage U = 1.5V and very small gap energy A — »• with 
the corresponding values for metallic SWNTs @, we 
obtain that the current values are around a factor six 
times smaller after the gap is opened. The crosses at 
the curves in Fig. 2 shown for A < hoj are calculated 
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FIG. 2: (Color online) Upper panel with the continuous 
curves shows the current-voltage curve obtained using ©-([91) 
for a nanotube of length L — 300 nm, l e — 1600 nm [2], |5|, 
and various gap energies A. The dashed curve is calculated 
for A — > and l e — 800nm. The crosses uses the full two 
electron band, three phonon type theory with exact transi- 
tion matrix elements. The dotted curves uses GW parameter 
values within the one-valley model for the electron-optical 
phonon scattering time [3(j. The lower panel shows the cor- 
responding momentum and position averaged backward scat- 
tering phonon distribution function n bs defined in Ref. 0. 



using the full two electron valley model coupled to Ei g 
T-phonons as well as A[ K-phonons where we used now 
the exact electron-phonon transition matrix elements in 
our numerics shown at the left hand side in (O and (J6]) 
for the r-phonons. We used here parameter values given 
in Refs. 0, M, EH and further that the phonon velocities 
of the transversal T-phonons is minus half of the velocity 
of the longitudinal ones [28| . 

In the lower panel in Fig. 2 we show the corresponding 
energy-position averaged backward phonon distribution 
function n bs ©-(HU). We obtain that the phonon 
distribution function is around three to four times smaller 
for a SWNT after gap opening at 1.5V when comparing 
it with the phonon distribution function values of the 
corresponding metallic nanotube [1, EH . 

Figure 3 shows the current- voltage characteristic of a 
metallic SWNT where just the gap is opened, i.e. A — > 0, 
for various nanotube lengths. Here we recorded more 
data points. We obtain a dip in the current- voltage char- 
acteristic with a large negative differential conductance 
for U > fiio/e by taking into account that hw ~ 0.1 7eV. 
Fig. 2 shows that this behavior survives even at much 
smaller elastic scattering lengths. We can understand 
this dip in the current curve in the following way: For 
small voltages ksT/e <C U <C hw/e, phonon scattering 
can be neglected and the total current is the sum of the 



FIG. 3: (Color online) We show the current-voltage curve 
for SWNTs with a gap A -> and certain SWNT lengths 
L (l e = 1600nm). The inset shows the current curve for a 
voltage U = IV and several nanotube lengths L as a function 
of the gap energy A. 



currents in each band separately. This current is for gen- 
eral A approximately given by 



denF^e) 



8e 



k B TUn 



A/2 



1 + exp 



A 



(13) 

with A = A/ksT and rip is the Fermi-function. At 
room temperature T « 300K and small gaps A —> 
we have Iq « 5.3/iA. We point out that Iq is indepen- 
dent of the length L of the SWNT. In ([13]) we neglect 
elastic scattering, which is only a small correction factor 
since L/l e < 1 for small SWNTs [H]. In the following 
discussion it is useful to have in mind that the energy 
broadening KbT of the electrons (holes) in the conduc- 
tion (valence) band is much smaller than the phonon 
frequency at room temperature, i.e. ksT/huj ~ 0.14. 
At higher voltages U > hw/e first interband scattering 
starts. Due to the electron distribution function factors 
in the phonon scattering terms in ([10)) and ([TTj) inter- 
band backward scattering from the conduction band to 
the valence band is then the dominant scattering contri- 
bution. This leads to a reduction of the current Iq of 
order ~ (e/h)kBTL/l ep where the electron-phonon 



mean free path is given by l ep = VFT e 



130nm hj 



At even higher voltages U ^ hw/e also innerband scat- 
tering starts. This leads to an energy diffusion of the 
electrons (holes) in the conduction (valence) band with 
an energy broadening from Ae ~ fc^T without scattering 
to Ae ~ (hw/7r)(l + L/l ep ) for eU > {hw/n){l + L/l ep ) 
[ll| . This broader energy band consists of equally spaced 
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subbands with separation distance Huj and energy width 
~ ksT. This follows from k B T <C hu and by taking 
into account the energy conservation in the electron- 
phonon scattering process. In this voltage regime the 
reduction of Iq due to interband scattering is changed 
tolf >w - (e/h)(k B T/n)(l + L/l ep )(hu/eU)L/l ep . This 
leads to a vanishing of interband scattering at high volt- 
ages shown in Fig. 3 as a dip in the current- voltage char- 
acteristic with negative differential conductance. From 
Fig. 3 we obtain that the considerations above are true 
only for L <C 500nm. For L w 500nm we obtain a satu- 
ration of the depth and the width of the voltage dip. The 
reason lies in the fact that at larger nanotube lengths also 
secondary scattering processes like interband scattering 
from the valence band to the conduction band has to be 
taken into account which was not done in the scaling con- 
siderations above. Also the phonon density dependence 
of the effective electron-phonon scattering length l ep 
should be included. 

At even higher voltages we obtain from Fig. 3 that 
lim.[/_>. 00 1 m Iq (apart from elastic scattering). This 
shows that a current reduction due to innerband scat- 
tering is very small in this regime. Phonons which 
are created by electron scattering in one band can only 
scatter electrons in the other band on a length Ax ~ 
L(l+ L / l cp )(hhj / ell) / n so that we can consider the Boltz- 
mann system (fT))- (fTTj) in each band separately at high 
voltages. This leads to the result that the electrons 
(holes) are effectively down (up) scattered in the con- 
duction (valence) band. Then, by taking into account 
the fact of existing band edges near the Fermi-niveau for 
small stresses, we obtain that innerband scattering does 
not influence the absolute current at high voltages. 

The inset in Fig. 3 shows the current- voltage curve for 
a SWNT at bias voltage U = 1 V for various nanotube 
lengths as a function of the gap energy A. The functional 
behavior of the reduction factor from the A — > current 
value is reproduced well by Eq. (fT3| . 

Recently, it was shown for graphene [2!| [3(| that 
screening effects in the electron-electron interaction 
are relevant for the calculation of the zone-boundary 
electron-phonon scattering time. In Ref. a numerical 
Green's function method (GW) was used in contrast to 
the standard density functional (DFT) theory in order to 
calculate this scattering time. Note that our one-valley 
model used so far is based on scattering parameters ob- 
tained by the latter calculation. The calculated values 
within the GW method are in better agreement with the 
experimentally obtained energy dispersion of the zone- 
boundary phonons for graphene than the values of the 
DFT calculation. For the zone-center phonons the agree- 
ment is less satisfactory with the experiments than in the 
DFT- method |3l|, |32(. To our knowledge up to now in 
contrast to the DFT calculations for SWNTs in Ref. 
an explicit GW calculation for the electron-phonon scat- 
tering time in carbon nanotubes is still missing. Since 
electron screening is highly dimension dependent this cal- 
culation would be of course of general interest. Neverthe- 



less we are now able to calculate from Ref. |3C§ within the 
GW method by using the zone-folding method [H effec- 
tive electron-phonon scattering times for zone-center and 
zone-boundary phonon scattering for SWNTs. By using 
Mathiessen rule we obtain an effective electron-phonon 
scattering length l cp = 68nm for an SWNT with diame- 
ter 2nm. Corresponding scattering times for the phonons 
due to electron-phonon scattering, the averaged phonon 
frequencies and phonon velocity for the one- valley model 
can be calculated accordingly by using the definitions 
in Ref. [iH . The acoustic-optical phonon scattering ra- 
tio T ac /pr op [l6j is then determined from the theoreti- 
cal determined current-voltage characteristic of the full 
two- valley model by comparison with the experimentally 
given current-voltage characteristic for the L — 300nm 
SWNT without stress. We obtain values T ac /pT op = 3.5 
for an optimal accordance to the differential conductance 
at high voltages. The absolute theoretical current- value 
is then around 20% lower in the high voltage regime than 
in the experiment. By using these theoretical GW quan- 
tities we show by the dotted curves in the upper panel in 
Fig. 1 the current- voltage characteristic for aL = 300nm 
SWNT and various gap energies A within an effective 
one- valley model. 

Up to now we only considered metallic nanotubes. 
Semiconducting nanotubes show already a gap A with- 
out stress. Since A 3> fhu for small diameter semicon- 
ducting SWNTs (~ 2nm) we obtain a vanishingly small 
current-voltage curve in this case. For semiconducting 
tubes with much larger diameters, the energy gap A is de- 
creasing but in order to obtain quantitatively the current- 
voltage characteristic one should take into account also 
the much larger electron-phonon scattering time r ep Q of 
such systems and further that now also higher transver- 
sal electron bands become important. Nevertheless, we 
expect from the discussion above a qualitative similar 
current- voltage behavior, with a negative differential con- 
ductance at U > huj/e as for the metallic nanotubes. 
Similar should be true for graphene nanoribbons. Re- 
member that the existence of a negative differential con- 
ductance region in the current-voltage characteristic is 
mainly based on the gap opening under stress and that 
we have k B T < fiw at room temperature. Here we 
note that corresponding to armchair nanotubes discussed 
above metallic zigzag-edged nanoribbons do not show a 
gap opening under axial stress in contrast to armchair 
edged ones. This is shown by numerical calculations in 
Ref. HH based on an atomistic model. This follows also 
immediately from the electron-phonon interaction Hamil- 
tonian for acoustic phonons in Ref. [TH which we used 
above in the derivation of |T]) and ([2|). 

Finally, we briefly discuss the current-voltage charac- 
teristic of a SWNT in the supcrclongation regime. In 
Ref. H we obtain that with increasing external stress de- 
fect stripes are created with increasing lengths for in- 
creasing stress. This is seen by a kink motion at the 
boundary surface of the SWNT in transmission electron 
microscope pictures Q. The number of these stripes is 
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proportional to the external stress whose lengths corre- 
sponds to missing lattice plane pieces in the crystal. In 
first approximation we expect that the current is given by 
the electrons running over complete lattice planes. The 
electrons which run over incomplete lattice planes are 
completely backscattered and these do not contribute to 
the current. For those electrons which run over the com- 
plete lattice planes we further have to consider the linear 
elasticity current corrections which we have discussed at 
length above. All this means that for armchair SWNTs, 
which do not show a stress dependent gap we expect a 
linear decreasing current curve as a function of the ex- 
ternal stress. Such a behavior is in fact experimentally 
observed in Ref . 

VI. CONCLUSION 

We have calculated the current-voltage characteristic 
of metallic carbon nanotubes under homogeneous ax- 
ial stress in the linear elastic regime lying on a sub- 
strate by means of a semi-empirical Boltzmann equa- 
tion for electrons and phonons. We have found that the 
stress leads to an energy gap which, causes a large neg- 
ative differential conductance at U > huj/e. Moreover 
this negative differential conductance is largest at small 



stress just after the gap is opened, getting smaller with 
a larger stress applied. Materials with such a behav- 
ior are very useful for the activation and deattenuation 
of oscillating circuits. For larger voltages at U ~ 1.5V 
and small stress just after the gap opening we obtain 
approximately six times smaller current values for the 
L = 300nm nanotube in comparison to the tube without 
stress. The corresponding momentum and position aver- 
aged optical phonon distribution function is around three 
to four times smaller. In addition to the discussion of the 
current-voltage characteristic in the linear elastic regime 
we have also discussed modifications of the characteristic 
in the superelongation regime. 

Concerning the current-voltage characteristic of large 
diameter semiconducting nanotubes and graphene 
nanoribbon, we have argumented that a qualitative simi- 
lar differential negative behaviour as for the metallic nan- 
otube should be seen. 
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